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INTRODUCTION 


Mie  theory  provides  a  means  of  calculating  the  absorp¬ 
tion  and  scattering  of  electromagnet j c  radiation  by  spheres 
of  any  size  compared  with  the  wavelength.  Although  exact, 
performing  full  Mie  theory  calculations  is  often  cumbersome: 
it  requires  an  excessive  amount  of  numerical  storage  space 
and  computer  time,  especially  when  the  particles  are  large 
compared  with  the  wavelength.  Under  certain  circumstances 
however,  it  is  possible  to  obtain  a  reasonable  estimate  of 
the  scattering  and  absorption  properties  of  spheres  without 
resorting  to  detailed  Mie  theory  calculations.  These  alter¬ 
native  approaches  to  exact  Mie  theory  calculations  are 
useful  provided  the  user  realizes  their  limitations. 

The  regions  of  validity  for  Mie  approximations  are 
usually  characterized  in  terms  of  the  size  parameter, 
X  =  2TTr/X,  where  r  is  the  radius  of  the  particle  and  \  is 
the  wavelength  of  the  impinging  radiation.  Unfortunately, 
the  use  of  an  approximation  over  Mie  theory  is  not  a 
clearcut  decision;  that  is,  the  accuracy  of  most 
approximations  depend  not  only  on  X,  but  to  some  extent  on 
the  complex  index  of  refraction  of  the  particle,  m  =  n  +  ik. 
Also  the  user  must  consider  the  degree  of  accuracy  required 
for  his  or  her  particular  problem. 

1.1  Organization  of  the  Report 

It  is  the  aim  of  this  report  to  explore  some  of  the  al¬ 
ternatives  to  Mie  theory  and  to  access  their  ranges  of 
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applicability.  Chapter  2  discusses  the  use  of  ray  optics  to 
obtain  the  phase  matrix  elements  of  a  sphere  that  is  large 
compared  with  the  wavelength.  Chapter  3  discusses  the  use 
of  complex  angular  momentum  theory  to  calculate  the 
efficiency  factors  for  extinction,  absorption  and  radiation 
pressure.  Chapter  4  compares  various  Mie  approximations  in 
the  scientific  literature  and  assesses  their  regions  of 
validity.  Chapter-  5  summarizes  the  main  results  of  this 
report.  For  reference.  Appendix  A  gives  the  equations  for 
the  various  Mie  approximations  investigated  in  this  report. 
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2.  THE  DSE  OF  RAY  OPTICS  TO  OBTAIN  THE  SCATTERING  MATRIX 
ELEMENTS  FOR  A  SPHERE 

2 . I  Overview 

The  angular  scattering  properties  of  spheres  can  be 
calculated  using  the  scheme  of  ray  optics  (van  de  Hulst^  and 
Liou  and  Hansen  )  provided  four  conditions  are  satisfied: 

1)  the  size  parameter  is  much  greater  than  one 

2)  the  real  part  of  the  index  of  refraction  is  greater 
than  one 

3)  the  user  is  interested  in  the  scattering  pattern  for 
a  small  spread  of  particle  sizes 

4)  the  far  field  scattering  pattern  is  desired. 

Tn  ray  optics,  incident  radiation  is  treated  as  local¬ 
ized  rays  whose  encounters  with  the  particle  are  governed  by 
the  laws  of  reflection  and  refraction  for  geometrical  op¬ 
tics.  The  localized  rays  can  be  reflected  externally  or  can 
enter  the  particle  where  they  are  either  absorbed  or 
internally  reflected  a  number  of  times  before  exiting. 
Also,  light  rays  passing  very  close  to  the  particle  are 
affected  by  Fraunhofer  diffraction.  Thus  the  total 

intensity  at  a  particular  scattering  angle  is  equal  to  the 
contribution  from  Fraunhofer  diffraction,  plus  that  for  rays 
emerging  from  the  particle  after  reflection  and  refraction. 

1.  van  de  Hulst,  H.  C.  (1957)  Light  Scattering  by  Small 

Particles ,  Wiley,  Inc.,  New  York,  470  pp. 

2.  Liou,  K.,  and  Hansen,  J.  E.  (1971)  Intensity  and  polar¬ 

ization  for  single  scattering  by  polydisperse  spheres: 
A  comparison  of  ray  optics  and  Mie  theory,  J.  Atmos. 

Sci .  28:995-1004. 
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(Strictly  speaking  it  is  not  correct  to  simply  add  the  in¬ 
tensities  of  the  separate  rays.  For  a  single  particle, 
emergent  rays  interfere  with  each  other  because  of  phase 
changes  caused  by  external  and  internal  reflection,  and  by 
differing  optical  paths.  These  interference  effects  lead  to 
rapid  oscillations  in  the  scattered  intensity  as  a  function 
of  scattering  angle  which  can  only  be  predicted  by  Mie  the¬ 
ory.  However  if  a  small  spread  of  particle  sizes  is  consid¬ 
ered,  any  phase  effects  are  averaged  out  when  integrating 
over  the  spread  of  particle  sizes.) 

In  ray  optics,  the  gain,  G,  is  used  to  describe  the 
relative  scattered  intensity  as  a  function  of  scattering  an¬ 
gle.  It  represents  the  ratio  of  scattered  intensity  to  that 
which  would  exist  if  it  were  scattered  isotropically  and 
conservatively^.  Thus  the  gain  is  defined  such  that 

1 

47r 

where  il  is  the  solid  angle  and  w^,  the  single  scattering 
albedo,  is  the  ratio  of  the  scattered  energy  to  that  scat¬ 
tered  and  absorbed  by  the  particles. 

Below  the  equations  governing  ray  optics  are  presented. 
These  equations  are  essentially  the  same  as  those  found  in 
van  de  Hulst^  and  in  Liou  and  Hansen^.  Next,  a  Fortran  sub¬ 
routine  called  RAYOPT  is  described  which  determines  the  an¬ 
gular  dependent  phase  matrix  elements  i^^,  i  2,  i^  and  i^. 
These  matrix  elements  relate  the  polarization  and  the  ampli¬ 
tude  of  the  incident  field  to  that  of  the  scattered  field. 


/(! 


U)d[! 
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The  results  from  RAYOPT  are  then  verified  against  those  from 
Liou  and  Hansen,  and  from  full  Mie  theory. 


2.2  Fraunhofer  Diffraction 

Diffraction  of  the  incident  beam  occurs  because  the 
wavefront  is  partially  obstructed  as  it  passes  /ery  close  to 
the  particle.  For  spheres,  the  gain  as  a  function  of  scat¬ 
tering  angle  resulting  from  diffraction  is  given  by 


fl,X)  =  G2^(  d,X) 


2X‘ 


(X  sin0) 
sin  9 


(2) 


where  the  superscript  F  denotes  the  gain  due  to  Fraunhofer 
diffraction,  is  the  first  order  Bessel  function  and  the 
subscripts  1  and  2  refer  to  the  gain  components  perpendicu¬ 
lar  and  parallel,  respectively,  to  a  particular  scattering 
plane.  The  diffracted  pattern  is  radially  symmetric  about 
the  forward  direction  and  becomes  more  concentrated  in  the 
forward  direction  as  the  size  parameter  increases.  For  the 
case  when  0=0,  Gj^^(0,X)  =  G2^(0,X)  =  0.5X^. 


2 . 3  Geonetrical  Optics 

The  geometry  of  the  problem  is  given  in  Figure  1.  The 
angle  between  an  incident  ray  and  the  tangent  of  the  local 
particle  surface  is  denoted  by  t  .  For  grazing  incidence, 
T  equals  0  degrees;  for  central  incidence,  r  equals  90  de¬ 
grees.  As  an  incident  ray  encounters  the  particle,  some  of 
its  energy  is  externally  reflected  while  the  rest  enters  the 
particle.  The  ratio  of  externally  reflected  energy  to  the 
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Figure  1.  Path  of  a  Light  Ray  According  to  Geometrical 
Optics,  Where  p  Is  the  Order  of  the  Ray  (see  Eq.  6) 

incident  energy  is  equal  to  the  squares  of  the  Fresnel  re¬ 
flection  coefficients,  r^  and  r2: 

sinr  -  m  sinr '  ^  ( 3  ) 

sinr  +  m  sinr* 

and 

2 

m  sinr  -  sinr'  (4) 

m  sinr  +  sinr' 

where  r'  is  the  direction  of  the  ray  upon  entering  the  par¬ 
ticle.  The  angle  r'  is  related  to  the  incident  angle,  t, 
through  Snell's  law 

cosr'  =  cosr/n  . 

Because  energy  is  conserved,  the  energy  entering  the  parti¬ 
cle  must  either  be  absorbed  or  emerge  after  any  number  of 
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internal  reflections.  For  a  nonabsorbing  particle,  the 

fraction  of  energy  emerging  after  refraction  with  no  inter- 

2  2 

nal  reflection  is  equal  to  (1-1  r^^  2 1  )  ?  after  one  internal 

2  2  2 

reflection,  it  is  I  r  2I  (l-lr^^  2^  ^  »  etc.  Here  it  is  con¬ 

venient  to  define 

2  ^^12^^  p  =  0 

^  r  ^  M  —  Ir-  I  ^  \  ^  I  I  -r-  |2\p— 1  £  .1  (6) 

where  p  =  0  means  external  reflection  and  p  >.  ^  means  re¬ 
fraction  plus  p-1  internal  reflections  (see  Fig.  1).  The 

2 

gain  for  geometrical  optics,  as  given  by  Liou  and  Hansen  , 
is  then 


(p) 


1,2^^'^^  -  2e  2^2^  exp(-4X  k  p  sinr'). 


(7) 


where  D  is  the  divergence,  defined  as 


SinrcoSr 


D  = 


sin0 


dd 


d  T 


and 


d  6  tanr 

-  =  2  -  2p  - 

dr  tanr' 


(8) 


(9) 


In  Eq.  7,  the  exponential  term  accounts  for  absorption  with¬ 
in  the  particle.  The  total  gain  at  a  particular  scattering 
angle  is  then 


G 


F 

1/2 


+ 


z 


p=0 


(10) 
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where  in  practice,  a  value  of  3  for  N  is  sufficient  to  ac¬ 
count  for  most  of  the  energy  scattered  by  geometrical  op¬ 
tics  . 


2.4  Relationship  Between  Gain  And  The  Phase  Matrix  Elements 

The  angular  dependent  phase  matrix  elements  relate  the 
polarization  and  amplitude  of  the  scattered  field  to  the  in- 
cident  field.  In  the  work  of  Blattner  ,  the  incident  and 
scattered  fields  are  described  by  the  parameters  I2,  U 
and  V  which  are  similar  to,  but  not  the  same  as,  the  Stokes 
parameters.  For  a  spherical  scatterer,  the  scattered  field 
(s  subscript)  is  related  to  the  incident  field  (o  subscript) 
through  the  amplitude  scattering  matrix 


'si 

i^  0  0  0 

'01 

's2 

X  2 

0  12  0  0 

'02 

'^s 

■  4.2r2 

0  0  i^  i^ 

«o 

V 

0  0  -i.  i. 

V 

s 

4  j  ' 

0 

where  R 
and  i  , 
matrix 

^1 
^2 
^3 
"4  ^ 


is  the  distance  from  the  scatterer  to  the  observer 
i2,  i3  and  i^  are  the  angular  dependent  scattering 
elements  defined  as 


=  clS^I  2 
=  cl  82!^ 

=  c  Re{Sj^S2*} 

=  -c  Im{S^S2*} 


3.  Blattner,  W.  (1972)  Utilization  instructions  for  opera¬ 
tion  of  the  Mie  programs  on  the  CDC-6600  computer  at 
AFCRL ,  F19628-70-C-0156,  Research  Note,  RRA-N7240,  Ra¬ 
diation  Research  Associates,  Inc.,  Fort  Worth,  Texas. 
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where  and  S2  are  the  familar  complex  amplitude  functions 

and  c  is  a  normalization  constant  such  that 


=  A 


{W)dJ2 


Coniparing  Equations  1  and  13  suggests  that 

11  =  G^/47rw^ 

12  = 

13  =  (Gj^G2)^/^/47rw^ 


2.5  Ray  Optics  Subroutine 

A  subroutine  called  RAYOPT  has  been  developed  to  deter¬ 
mine  the  phase  matrix  elements  i^,  i2,  i3  and  i^  at  user 
specified  scattering  angles.  The  inputs  to  the  subroutine 
are  listed  in  Table  1.  A  flow  diagram  for  RAYOPT  is  given 
in  Figure  2  and  a  brief  description  of  the  important  vari¬ 
ables  and  functions  is  given  in  Table  2. 

For  geometrical  optics,  RAYOPT  determines  the  gains  at 
the  user  scattering  angles  in  the  following  manner.  With  p 
held  fixed,  the  total  deviation  of  a  ray  from  its  original 
path,  6',  is  computed  for  incident  angles,  7,  between  0  and 


90  degrees  in  steps  of  0.5  degrees  where 


-  2  T  -  2pT' 


(15) 


Next,  each  user  scattering  angle  is  transformed  so  that  it 
also  represents  the  total  deviation  of  a  ray  from  its  origi¬ 
nal  path.  This  deviation  is  then  compared  with  the  devia- 
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Table  1 .  Input  Parameters  for  RAYOPT 


1 

VARIABLE 

DESCRIPTION 

NP 

Maximum  Value  for  p  When  Performing  the 

Ray  Tracing 

NTHETA 

Number  Of  User  Scattering  Angles.  The 

Maximum  is  150 

PM 

Complex  Index  of  Refraction  of  the  Particle. 
The  Imaginary  Part  Is  Positive 

THETdSO) 

Array  Containing  the  User  Scattering  Angles 

WO 

Single  Scattering  Albedo  of  the  Particle 

X 

Size  Parameter 

tions  for  known  incident  angles  to  determine  what  incident 
angles  yield  emergent  rays  at  a  user  scattering  angle.  When 
a  user  scattering  angle  falls  between  two  values  of  0 ' ,  the 
incident  angle  is  obtained  by  means  of  a  linear  interpola¬ 
tion.  With  the  incident  angle  now  known,  the  gain  at  the 
user  scattering  angle  is  calculated  using  Eqs .  3-9. 

Before  computing  the  gain,  RAYOPT  checks  to  make  sure 
the  incident  angle  does  not  lead  to  either  a  rainbow,  i.e. 
d6/dT*~0,  or  a  glory,  i.e.  6=  180  and  sin  6  =  0.  (The  scheme 
of  ray  optics  predicts  infinite  intensities  at  these  angles, 
see  Eq.  8.)  If  a  rainbow  or  glory  angle  is  found,  the  gain 
calculation  is  performed  for  an  incident  angle  given 
by  T =  r -  2  degrees.  The  resulting  gain  is  then  used  unless 
it  is  still  excessively  large  OlOO);  in  that  case,  the  gain 
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Figure  2 


Comput*  th«  Qain 
at  tha  User  Scattaring  Anglaa 
Du#  to  Difiractlon 

i 

Calculata  tha  QaIn 
at  tha  User  Scattarlng  Angle# 

Dua  to  External  Raflactlon  (P  ■  0) 

i 

Oatarmlne  the  Qain 
at  tha  Uear  Scattaring  Anglea 
lor  All  P  That  Do  Not  Have 
a  Rainbow  or  Qlory  (Criterion:  P  <  n') 

i 

Compute  tha  Qain 
at  tha  User  Scattering  Anglee 
tor  All  Remaining  P  Upto  NP 

i 

Convert  tha  Qalna  to  Matrix  Eiemanta 

i 

Return 


Flow  Diagram  of  the  Subroutine  RAYC'PT 
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Table  2.  Important  Variables  and  Functions  in  RAYOPT 


PARAMETER 

DESCRIPTION 

GAINP 

Function  That  Returns  the  Gain  for  a 
Given  Incident  Angle  and  p 

GAINKISO)  , 
GAIM2(150) 

Keeps  a  Running  Total  of  the  Gains 
at  the  User  Scattering  Angles 

GDIPPdSO) 

Gains  at  the  User  Scattering  Angles 

Due  To  Fraunhofer  Diffraction 

J1 

Function  That  Returns  The  Bessel 
Function  Of  The  First  Kind 

PNAUT 

Maximum  Value  of  p  that  Does  Not  Have 
a  Rainbow  or  Glory 

RENTl (150) , 
RENT2(150) 
ENT3(150) , 
ENT4(150) 

Angular  Dependent  Matrix  Elements 
^2,  ^4  Respectively 

TAG 

Angle  of  An  Incident  Ray  With  Respect 
To  The  Particle  Surface 

THETPRdSl) 

Array  Containing  the  Total  Deviation 
of  Emergent  Rays  for  Incident  Rays 
between  0  and  90  degrees  in  Steps  of 

0 . 5  degrees 

12 


is  set  equal  to  100.  The  gain  calculations  are  simplified 
for  p  _<  n  because  they  do  not  contain  rainbows  or  glorys. 
These  values  of  p  are  done  first  in  RAYOPT  where  the  amount 
of  computer  time  is  reduced  significantly.  Finally  for  cen¬ 
tral  incidence,  where  t  is  90  degrees  and  0  is  0  degrees, 
the  gain  equals  2  =  0.5(e^  (l-p/rx)^  )exp{-AX  k  p)  . 

2.6  Saaiple  Calculations  Using  RAYOPT 

To  ensure  that  RAYOPT  is  working  properly,  a  series  of 
calculations  have  been  performed  and  compared  with  those 
from  other  researchers.  Figures  3  and  4  duplicate  Figures  1 


Scattering  Angle  (6)  Scottering  Angle  {&) 


2 

Figure  3.  Recalculation  of  Figure  1  of  Liou  and  Hansen  , 
Angular  Gain  for  Geometrical  Optics  using  RAYOPT.  The  Ver¬ 
tical  Scale  Applies  to  the  Lowermost  Curves  (n  =  1.45  and 
2.00)  while  the  Other  Curves  Are  Displaced  Upwards  by  Fac¬ 
tors  of  10^ 
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Scattering  Angle  (0)  Scottering  Angle  (0) 

2 

Figure  4.  Recalculation  of  Figure  2  of  Liou  and  Hansen  , 
Degree  of  Polarization  for  Geometrical  Optics  using  RAYOPT. 
The  Vertical  Scale  Applies  to  the  Lowermost  Curves  (n  =  1.45 
and  2.00)  while  the  Other  Curves  Are  Displaced  Upwards  by 


2 

and  2  of  Liou  and  Hansen  .  There  is  complete  agreement  ex¬ 
cept  in  some  instances  at  the  rainbow  and  glory  angles:  Liou 
and  Hansen  did  not  say  how  they  dealt  with  these  angles. 

Figures  5  and  6  give  comparisons  between  ray  optics  and 
Mie  theory  for  the  matrix  elements  i^^,  i2  and  i^.  Figure  5 
represents  the  scattering  by  a  particle  with  an  index  of  re¬ 
fraction  of  1.33  +  O.Oi  while  in  Figure  6,  an  index  of  re¬ 
fraction  of  1.'^  +  O.Oli  has  been  used.  In  both  figures, 
the  matrix  elements  for  ray  optics  represent  a  size  parame¬ 
ter  of  400.  When  performing  the  Mie  calculations,  it  was 
necessary  to  use  a  particle  number  density  distribution, 
N(r),  described  by  a  log-normal  distribution  because  inter¬ 
ference  effects  lead  to  rapid  oscillations  in  the  angular 
scattering  pattern  of  a  single  sphere.  This  analytical  for¬ 
mulation  is  given  by 


dN  { r ) 
dlog ( r ) 


(27r)l/2loga 


exp  - 


log (r/r^) 
2  log  a 


(16) 


where  r^,  a  and  are  the  geometric  mean  radius,  geometric 
standard  deviation  and  total  number  of  particles  respec¬ 
tively.  In  Figures  5  and  6,  an  r^  of  63.66  micro¬ 
meters,  a  of  0.1,  and  of  1  particle  cm~^  have  been  used. 
The  value  for  r^  was  chosen  because  it  converts  to  a  size 
parameter  of  400  for  1  micrometer  radiation.  The  range  of 
the  Mie  calculations  was  from  50  to  100  micrometers. 

Figure  5(a)  shows  good  agreement  between  the  values  of 
for  ray  optics  and  Mie  theory  except  for  the  rainbow  re¬ 
gion  from  130  to  140  degrees  and  at  180  degrees.  The  peak 
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Fiqure  5.  Values  of  (a)  i, ,  (b)  for  Ray  Optics 
and  Mie  Theory  Where  There^Is  No  Absorption  Within 

the  Particle 


m  =  1.54  +  O.Oli  - 

Roy  Optics 

X  =  400 

.  Mie  Theory 

1 

- -  Rotio 

1 

1 

V' 

JU  60  90  120  150  180 

Scattering  Angle  (0) 


Figure  6.  (continued)  Values  of  (c)  for  Ray 
Optics  and  Mie  Theory  Where  There  Is  Absorption 
Within  the  Particle.  In  (c)  the  Values  Beyond 
About  70  Degrees  Are  Not  Plotted  Because  i^  for 
Mie  Theory  Is  Negative 


at  180  degrees  for  Mie  theory  is  related  to  the  familiar 
glory  phenomenon.  Ray  optics  does  not  predict  the  glory  and 
here  lies  one  of  its  major  weaknesses.  Figure  5(b)  suggests 
only  fair  agreement  between  the  values  of  12  for  ray  optics 
and  Mie  theory.  In  Figure  5(c),  the  values  of  i^  for  ray 
optics  are  close  to  those  for  Mie  theory  provided  the  scat¬ 
tering  angle  is  less  than  80  degrees.  Beyond  this  angle 
however,  ray  optics  cannot  be  used  in  place  of  Mie  theory 
because  i^  for  Mie  theory  is  negative.  This  limitation  of 
ray  optics  is  not  as  severe  as  it  may  appear  since  li^l  is 
small  beyond  80  degrees. 

Figures  6(a)  and  6(b)  show  better  general  agreement  be¬ 
tween  Mie  theory  and  ray  optics  except  near  the  diffraction 
peak.  The  improved  accuracy  of  ray  optics  can  be  attributed 
to  the  absorption  within  the  particle  which  eliminates  the 
less  accurate  contributions  from  p  ^  1.  However  discrepan¬ 
cies  still  exist  in  the  values  of  ±2  for  scattering  angles 
between  60  and  80  degrees.  Here,  ray  optics  predicts  values 
for  £2  that  are  10  to  100  times  greater  than  Mie  theory.  In 
Figure  6(c),  the  values  of  i^  for  ray  optics  are  close  to 
those  for  Mie  theory  provided  the  scattering  angle  is  less 
than  60  degrees.  Beyond  60  degrees,  Mie  theory  predicts 
negative  values  for  i^. 
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3. 


CALCDLATING  NIE  EFFICIENCY  FACTORS  WITH  COMPLEX  ANGULAR 
MOMENTUM  THEORY 


3.1  Overview 

Three  parameters  often  obtained  from  Mie  theory  calcu¬ 
lations  are  the  efficiency  factors  for  extinction,  absorp¬ 
tion  and  radiation  pressure  which  are  refered  to  as 
Qabs  Qpr  ^respectively.  and  are  related  to  the 

extinction  and  absorption  cross  sections,  and  C^j^g,  by 

C. 


'ext 


■'ext 


(17) 


and 


'abs 


■abs 


(18) 


where  A  is  the  cross-sectional  area  of  the  sphere  and  the 
scattering  efficiency  follows  as  Qgj,^  “  ®ext  ~  ®abs’ 
radiation  pressure  efficiency,  which  measures  the  total  rate 
of  momentum  transfer  to  the  particle,  can  be  converted  to 
the  asymmetry  parameter,  g,  by 

^  ^  ^^ext  ~  ^pr^/^sct  • 


Complex  angular  momentum  theory  can  be  used  to  deter¬ 
mine  the  Mie  efficiency  factors  of  a  sphere  that  is  large 
compared  with  the  wavelength.  The  governing  equations  have 

4 

been  developed  by  Nussenzveig  and  Wiscombe  .  (These  equa¬ 
tions  have  been  replicated  in  Appendix  A  due  to  their  ex- 


4.  Nussenzveig,  H.  M.,  and  Wiscombe,  W.  J.  (1980)  Efficiency 
factors  in  Mie  scattering,  Phys.  Rev.  Letters,  45 ; 
1490-1494  . 
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treme  length.)  Nussenzveig  and  Wiscombe  have  shown  that 
when  compared  with  exact  Mie  values,  complex  angular  momen¬ 
tum  theory  gives  relative  errors  of  _(1-10)%  for  X  -  10  and 
-(0.01-0.001)%  for  X  =  1000.  The  theory  does  not  predict 
the  rapid  changes  in  the  values  of  and  as  a  func¬ 
tion  of  size  parameter,  but  this  high-frequency  "ripple 
structure"  can  be  ignored  in  many  applications. 

3 . 2  Software 

Source  code  has  been  obtained  from  Warren  Wiscombe^ 
which  calculates  Qgbs  ^pr  using  complex  angular 
momentum  theory.  The  code  was  written  in  standard  For¬ 
tran  77  and  therefore,  required  very  few  modifications  to 
operate  on  the  AFGL  CYBER  system.  Specifically,  to  avoid 
underflow  errors,  checks  have  been  inserted  before  calls  to 
the  Fortran  CEXP  function.  In  addition  the  ACCUR  parameter 
in  the  code,  which  establishes  a  convergence  criterion  for 
numerical  integrations,  has  been  changed  from  10  to  10  . 
This  change  was  recommended  by  Warren  Wiscombe  and  acts  to 
speed  up  the  program  without  a  significant  loss  of  accuracy. 

3.3  Timing  Comparisons  with  Mie  Theory 

It  is  well  known  that  the  computation  time  for  Mie  cal¬ 
culations  is  directly  proportional  to  the  size  parameter. 
Figure  7  illustrates  this  clearly  where  the  contours  repre¬ 
sent  CPU  times  for  various  size  parameters  and  imaginary  in¬ 
dices  of  refraction.  The  CPU  times  were  determined  using 

5.  Wiscombe,  W.  J.,  private  communication. 
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the  Fortran  SECOND  function. 


available  on  the  AFGL  Cyper 


computer . 

For  complex  angular  momentum  theory,  the  expression  for 
Qext  is  a  "one  line"  expression  that  requires  very  little 
computer  time.  On  the  otherhand,  the  expressions  for  Q^j^g 
and  involve  numerical  integrations  that  require  more  CPU 

time.  Therefore,  a  series  of  timing  runs  have  been  per¬ 
formed  for  a  wide  range  of  size  parameters  and  indices  of 
refraction  to  determine  where  complex  angular  momentum  the¬ 
ory  is  faster  than  Mie  theory.  In  the  comparisons,  the  CPU 
time  for  complex  angular  momentum  theory  represents  the  to- 

(Q^ 


tal  time  required  to  compute  Q 
-pr 


ext'  '^abs  ^pr*  '^ext' 

Qabs  Qpr  computed  simultaneously  during  a  single  Mie 

run.  ) 


Figure  8  gives  CPU  times  for  complex  angular  momentum 
theory  as  a  function  of  size  parameter  and  imaginary  index 
of  refraction.  The  figure  suggests  that  there  is  broad  min¬ 
ima  in  the  CPU  time  for  imaginary  indices  of  "(0.01-1). 
Other  timing  runs  indicated  that  the  CPU  times  were  not  in¬ 
fluenced  by  the  real  index  of  refraction. 

Figure  9  compares  the  CPU  times  for  complex  angular  mo¬ 
mentum  theory  against  those  for  Mie  theory  where  the  compar¬ 
ison  is  expressed  as 
T 

,  APX 

log  -  .  (20) 

T 

MIE 

In  Figure  9,  it  can  be  seen  that  complex  angular  momentum 
theory  is  faster  than  Mie  theory  for  size  parameters  greater 
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Innaginory  Index 


Size  Poronoeter 


Figure  8.  Computation  Time  for  Complex  Angular  Momentun 
Theory  (in  CPU  seconds)  as  a  Function  of  Size  Parameter  and 
the  Imaginary  Part  of  the  Index  of  Refraction.  The  Shaded 
Region  Denotes  Where  the  Computation  Time  Is  Less  Than  0.03 
Seconds 
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0.8 


Size  Parameter 


Figure  9.  Comparison  of  the  Computation  Times  for  Complex 
Angular  Momentum  Theory  Against  Those  for  Mie  Theory.  The 
Shaded  Region  Indicates  Where  Complex  Angular  Momentum 
Theory  Is  Faster  Than  Mie  Theory 


than  200  to  500,  depending  on  the  imaginary  index  of  refrac¬ 
tion.  Furthermore,  percent  errors  are  generally  less  than 
1%  for  size  parameters  greater  than  500  so  Mie  theory  can  be 
replaced  by  complex  angular  momentum  theory  without  sacri¬ 
ficing  accuracy. 

3.4  A  Singularity  in  the  Expression  for 

In  the  course  of  testing  the  source  code  provided  by 
Warren  Wiscombe,  it  was  discovered  that  a  singularity  exists 
for  whenever  the  imaginary  index  was  less  than  about 
10“^  and  the  real  index  equals  3,  5,  7...  Specifically, 
this  ill-behavior  arises  from  a  division  by  zero  in  the  sum¬ 
mation  term  of  Eq.  1  of  Nussenzveig  and  Wiscombe^.  The  ex¬ 
tent  of  the  ill-behavior  can  be  seen  in  Figures  10(a)  and 
10(b)  where  percent  errors,  when  compared  with  Mie  theory 
values,  are  given  for  real  parts  of  2.5  and  3.0.  Subsequent 
testing  showed  that  the  expression  for  can  be  used  con¬ 
fidently  provided 

I n  -  n* I  >  0 . 05  ( 21 ) 
where  n*  =  3 ,  5 ,  7 . .  . 


27 
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Figure  10.  Percent  Errors  in 
Momentum  Theory  for  (a)  n  =  2.50 


for  Complex  Angular 


(b)  10 


Figure  10.  (continued)  Percent  Errors  in  for  Complex 

Angular  Momentum  Theory  for  (b)  n  =  3.00.  In  (b) ,  the 
Shaded  Region  Illustrates  the  Singularity 
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4.  OTHER  APPROXIMATIONS  FOR  AND 

This  chapter  examines  various  approximations  for 

and  Q  ,  that  are  available  in  the  scientific  literature, 
'^abs 

Specifically,  regions  of  validity  are  established  for  each 
approximation  by  comparing  exact  Mie  value  against  the  cor¬ 
responding  approximate  value.  A  wide  range  of  size  parame¬ 
ters  and  indices  of  refraction  were  considered  in  the  analy¬ 
sis.  In  addition,  the  regions  of  validity  are  represented 
in  a  way  that  that  permits  the  various  approximations  to  be 
intercompared .  The  purpose  of  this  latter  test  is  to  see 
where  better  approximations  are  needed. 


4 . 1  Approximations  That  Were  Investigated 

Five  approximations  for  have  been  investigated  in 

c  n 

this  report:  two  for  small  particles  '  and  three  for  large 

4  8  9  .  . 

particles  '  '  .  For  a  total  of  seven  approximations 

were  looked  at:  those  given  by  plus  two  other  approx¬ 

imations  for  large  particles^® ' .  A  brief  description  of 


6.  Rayleigh,  Lord  (1871)  On  the  light  from  the  sky,  its  po¬ 

larization  and  colour,  Philos.  Mag. ,  £1:107-120,  274- 
279  (reprinted  in  Scientific  Papers  by  Lord  Rayleigh, 
Vol.  1:1869-1881,  No.  8,  Dover,  New  York,  1964). 

7.  Wiscombe,  W.  (1980)  Improved  Mie  scattering  algorithms, 

Appl .  Opt . ,  £:1505-1509. 

8.  Deirmend jian,  D.  (1969)  Electromagnetic  Scattering  on 

Spherical  Polydispersions,  Elsevier,  New  York. 

9.  Ackerman,  S.  A.  and  Stephens,  G.  L.  (1987)  The  absorption 

of  solar  radiation  by  cloud  droplets:  an  application 
of  anomalous  diffraction  theory,  J .  Atmos .  Sci . , 
£4:1574-1688. 

10.  Bohren,  C.  F.  and  Nevitt,  T.  (1983)  Absorption  by  a 

sphere:  A  simple  approximation,  Appl .Opt . ,  22 : 774-775  . 
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how  each  approximation  works  is  given  in  Table  3.  For  ref¬ 
erence,  the  equations  governing  the  various  approximations 
have  been  reproduced  in  Appendix  A  of  this  report. 

4 . 2  Conparison  Scheme 

4.2.1  Realm  of  Comparison 

Table  4  summarizes  the  size  parameters  and  complex  in¬ 
dices  of  refraction  that  are  investigated  in  the  current 
analysis.  These  values  encompass  the  realm  of  particles 
normally  encountered  in  atmospheric  aerosol  modeling.  For 
each  approximation,  n  is  held  fixed  and  the  exact  Mie  and 
approximate  values  are  calculated  at  incremental  steps  of 
the  size  parameter,  X,  and  the  imaginary  index  of  refrac¬ 
tion,  k.  The  steps  are  equally  spaced  in  the  logarithms  of 
X  and  k  where  the  spacing  consists  of  fifty  steps  per  de¬ 
cade  of  X  and  twenty  five  steps  per  decade  of  k.  Similar 
sets  of  data  are  also  computed  for  other  fixed  values  of  n. 

4.2.2  Smoothing  Function 

Before  constructing  error  plots,  data  for  Mie  theory 
and  the  approximations  have  been  convolved  with  an  eleven 
point  smoothing  function,  S.  Specifically,  the  smoothing 
function  is  passed  over  the  data  for  X  1  while  k  is  held 


11.  Levine,  P.  H.  (1978)  Absorption  efficiency  for  large 

spherical  particles;  A  new  approximation,  Appl .  Opt . , 
17:3861-3862. 
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Table  3.  Summary  of  the  Approximations  That  Were  Investigated 
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Table  Particle  Characteristics  That  Were 
Investigated 


PARAMETER  RANGE 


Real  Part  of  the 

1.1  -  ; 

3.0 

Index  of  Refraction,  n 

Imaginary  Part  of  the 

10-6  . 

10 

Index  of  Refraction,  k 

Size  Parameter,  X 

10“^  - 

10 

fixed  so  that 
10 


.-i: 

L=0 


(22) 


i+L-5 


where  represents  the  Mie  or  approximate  value  at  conse¬ 
cutive  steps  of  the  size  parameter,  X^ .  The  weighting  fac¬ 
tors  for  the  smoothing  function  are  given  by 


(“) 


102^ 


,  L  —  0,  10 


(23) 


where 


is  the  binomial  coefficient, 

10!/((L! )(10-L) ! )  . 


(24) 


The  smoothing  function  acts  to  eliminate  the  high-fre¬ 
quency  ripple  structure  that  occurs  in  the  exact  Mie  values 
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for  the  efficiency  factors  as  a  function  of  X  for  fixed  k. 
For  reference,  the  ripple  structure  in  is  illustrated 

in  Figure  11  for  fixed  values  of  k.  The  reason  for  elimi¬ 
nating  the  ripple  structure  is  twofold: 

1)  it  allows  for  a  fairer  assessment  of  the  various  ap¬ 
proximations  because  the  ripple  structure  can  only 
be  calculated  with  Mie  theory 

2)  since  many  aerosol  modeling  applications  involve  a 
polydispersion  of  particles,  the  effects  of  the  rip¬ 
ples  are  "washed  out"  when  integrating  over  particle 
size . 

The  width  of  the  smoothing  function  has  been  chosen  to  pre¬ 
serve  the  large  scale  oscillations,  known  as  interference 
structure . 


4.2.3  Percent  Error  Contour  Plots 

For  each  approximation,  two-dimensional  percent  error 
diagrams  have  been  developed  for  fixed  values  of  n.  The 
percent  errors  are  computed  by 


^X.,k.  ^^®^^X.,k.  ^X.  ,k  ,k  . 

I  D  D  1-  D  I  D 


(25) 


where  A  and  M  are  the  (smoothed)  approximate  and  Mie  values, 
respectively,  evaluated  at  size  parameter  X^  and  imaginary 
index  of  refraction  kj.  Contour  lines  have  been  drawn  for 
percent  errors  of  8,  20  and  50%. 


4.3  Results 


4.3.1  Regions  of  Validity  for  the  Extinction  Approximations 

Figures  12(a)  -  12(e)  give  representative  percent  error 
diagrams  for  each  of  the  extinction  approximations.  Each 
plot  is  for  n  =  1.54,  and  the  regions  where  percent  errors 
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Figure  12,  (continued)  Percent  Errors  in  for 

Approximation  of  (b)  Wiscombe,  for  n  =  1.54 
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Size  Parameter 


Figure  12,  (continued)  Percent  Errors  in  Q  for  the 
Approximation  of  (c)  Deirmend jian,  for  n  = 


Size  Parameter 


Figure  12.  (continued)  Percent  Errors  in  Q  for  the 
Approximation  of  (d)  Ackerman  and  Stephens,  for  n  =  1.54 
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Figure  12.  (continued)  Percent  Errors  in  Q  for  the 
Approximation  of  (e)  Nussenzveig  and  Wiscombi,  for  n  =  1.54 
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are  less  than  8%  have  been  shaded.  For  brevity,  plots  for 
other  values  of  n  have  not  been  included  here. 


Figures  12(a)  -  12(e),  and  similar  plots  for  other  n, 
have  been  used  to  determine  the  regions  of  validity  of  the 
various  approximations  for  These  results  are  given  in 
Figure  13.  The  dashed  line  for  the  Nussenzveig  and  Wiscombe 
approximation  indicates  the  useful  range  if  the  singularity 
that  arises  for  odd  integer  values  of  n  is  avoided  (see  Sec¬ 
tion  3.4). 


In  general,  existing  approximations  for  do  reason¬ 
ably  well  for  both  large  and  small  size  parameters.  Of  the 
large  particle  approximations,  the  Ackerman  and  Stephens  ap¬ 
proximation  extends  to  the  lowest  size  parameter.  If  higher 
accuracies  are  needed  (i.e.  0.1  to  1%  or  better),  the 
Nussenzveig  and  Wiscombe  approximation  has  a  wider  useful 
range  than  the  other  large  particle  approximations  because 
it  reproduces  the  higher-order  ripple  structure  better. 
Most  approximations  however,  have  trouble  predicting  the 
magnitude  and  location  of  the  first  interference  peak  which 
occurs  between  size  parameters  of  1  and  10  (see  Figure  11). 
On  the  other  hand,  percent  errors  within  this  narrow  band 
are  generally  less  than  20%  for  many  of  the  approximations. 


4.3.2  Regions  of  Validity  for  the  Absorption  Approximations 

Figures  14(a)  -  14(g)  give  representative  percent  error 
diagrams  for  each  of  the  absorption  approximations.  Each 
plot  is  for  n  =1.54,  and  the  regions  where  percent  errors 
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Figure  13.  Regions  of  Validity  for  the  Mie  Approximations 
for  The  Shaded  Regions  Indicates  Where  All  of  the 
Approximations  Give  Percent  Errors  >  8%.  The  Dashed  Line 
for  the  Nussenzveig  and  Wiscombe  Approximation  Indicates  the 
Useful  Range  If  the  Singularity  That  Arises  for  Odd  Integer 
Values  of  n  Is  Avoided 
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Figure  14.  Percent  Errors  in  Q  , 
(a)  Rayleigh,  for  n  =  1.54 


for  the  Approximation  of 
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Figure  14.  (continued)  Percent  Errors  in  Q  for  the 
Approximation  of  (b)  Wiscombe,  for  n  =  1.54 
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Figure  14.  (continued)  Percent  Errors  in  for  the 

Approximation  of  (d)  Ackerman  and  Stephens,  for  n  =  ] . 54 
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Figure  14.  (continued)  Percent  Errors  in  Qgj^g  for  the 
Approximation  of  (g)  Levine,  for  n  =  1.54 
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are  less  than  20%  have  been  shaded.  For  brevity,  plots  for 
other  n  have  not  been  included  here. 

Figures  14(a)  -  14(g),  and  similar  plots  for  other  n, 
have  been  used  to  determine  the  regions  of  validity  of  the 
various  approximations  for  Q^j^g  •  These  results  are  given  in 
Figure  15.  The  region  of  validity  of  the  Nussenzveig  and 
Wiscombe  approximation  is  drawn  with  a  dashed  line  because 
it  is  only  computionally  efficient  for  size  parameters 
greater  than  about  200,  relative  to  exact  Mie  theory  (see 
Section  3.3),  where  it  is  accurate  to  about  0.1%.  Also,  to 
avoid  clutter,  the  region  of  validity  of  the  Aclcerman  and 
Stephens  approximation  has  not  been  shown;  however,  its 
omission  does  not  alter  the  shape  of  trie  shaded  region  in 
Figure  15  because  the  Aclcerman  and  Stephens  approximation 
for  is  similar  to  the  Bohren  and  Nevitt  approximation. 

Figure  15  suggests  that  existing  approximations  for 
Qabs  less  accurate  than  those  for  For  example, 
the  tolerable  error  has  been  relaxed  to  20%  because  signifi¬ 


cant  errors  extend  to  rather  large  size  parameters.  Of  the 
large  particle  approximations,  the  Bohren  and  Nevitt  approx¬ 
imation  appears  to  have  the  largest  region  of  applicability. 
In  contrast  to  Figure  13,  the  percent  errors  within  the 
shaded  region  of  Figure  15  often  exceed  50%. 
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Figure  15.  Regions  of  Validity  for  the  Mie  Approximations 
for  Qgjjg*  The  Shaded  Regions  Indicates  Where  All  of  the 
Approximations  Give  Percent  Errors  >  20%.  The  Ackerman  and 
Stephens  Approximation  Is  Omitted  for  Clarity,  But  Is 
Similar  to  That  of  Bohren  and  Nevitt.  The  Nussenzveig  and 
Wiscombe  Approximation  Is  Shown  with  a  Dashed  Line  Because 
It  Is  Only  Computationally  Efficient  for  X  >  200 
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5.  SUMMARY  AND  CONCLUSIONS 


This  report  has  explored  alternative  algorithms  to  Mie 
theory.  Specifically,  Chapter  2  discussed  the  use  of  ray 
optics  to  calculate  the  phase  matrix  elements  of  polydis- 
perse  spheres  in  geometrical  optics  limit.  Generally  speak¬ 
ing,  ray  optics  should  be  used  for  size  parameters  greater 
than  400,  and  becomes  more  accurate  for  particles  containing 
significant  absorption.  Chapter  3  explored  the  use  of  com¬ 
plex  angular  momentum  theory  to  determine  Qgbs  ^pr 
(or  g) .  It  was  shown  that  although  accurate  for  size  param¬ 
eters  greater  than  20,  this  technique  is  only  computation¬ 
ally  efficient  for  size  parameters  greater  than  500.  In 
Chapter  4,  approximations  for  and  the  scien¬ 
tific  literature  were  investigated  and  their  regions  of  va¬ 
lidity  were  established.  It  was  shown  that  existing  approx¬ 
imations  can  give  reasonable  estimates  of  for  most  size 
parameters  and  complex  indices  of  refraction.  On  the  other 
hand,  the  approximations  for  are  much  less  accurate. 
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Appendix  A 


Equations  for  the  Various  Mie  Approxinatiuns 

Investigated  in  this  Report 


This  appendix  gives  the  equations  used  to  calculate  the 
extinction  efficiency,  Qg^t'  absorption  efficiency, 

Qabs'  approximations  discussed  in  this  report. 

For  some  of  the  approximations,  the  scattering  efficiency, 
Qsca»  calculated  in  place  of  either  or  Q^j^g  where 

the  three  terms  are  then  related  by  =  Q  __  +  Q  k  •  The 

list  of  equations  serves  only  as  an  overview  and  the  reader 
should  consult  the  original  references  for  more  detailed 
information.  Also,  some  of  notation  from  the  original  ref¬ 
erences  has  been  changed  such  that  it  is  consistent  with  the 
terminology  used  in  this  report.  In  the  equations  below,  X 
represents  the  particle  size  parameter  and  m  =  n  +  ik  repre¬ 
sents  the  complex  index  of  refraction  of  the  particle. 
Equations  for  the  Rayleigh  Approxiaation^ 

The  Rayleigh  approximation  for  Qgj^g  is  computed  as 


m-'-  1 


m‘  +  2 


and  Q_„.  is  given  by 

SCct 


o  =  — 


m 


1 


+  2 


where  -  Qg^a  ^abs 


(A-1) 


(A-2) 


A-1 
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lations  for  the  Wiscombe  Approximation 


The  Wiscombe  approximation  for  Q  ^  is  computed  as 


0...  =  + 


(A-3) 


and  is  given  by 


Q  =  6X  T 

'^sca 


{A-4) 


Qabs  =  Qext  -  Qsca 


1  y' 4. 

-  ,  .  -1  ^~10  ^  1400  ^ 

a  1  —  2 1  3 


(A-5) 


.  2  2  ,1  + 

1  _2m^-5  2 

^  30  ^ 


(A-6) 


ix^ 

15 


2m^  +  3  -  — -  ~  ^X^ 
14 


(A-7) 


D  =  m'+2+  (l  -  x"-  ■^— 350  ^ 

V  10  /  1400  ^ 


+  2i 


■  m^-1  „3 


3  /  12' 


(A-8) 


.  i2  A  2  c  , 


(A-9) 


Ar2 


Q 

Equations  for  the  Deinnendiian  Approximation 

The  Deirmend jian  approximation  for  is  computed  as 


^ext  ^VHext^^ 

(A-10) 

and  Qgt)s  given  by 

Q  =  Q  K  (1  +  d) 

'^abs  '^VHabs 

(A-11) 

where  Qyjjext  extinction  efficiency 

Hulst^, 

from 

van  de 

4  cos  P 

^vHext-^  p  exp(-p  tan  p)  sin(p-p) 

2 

/  QOS 

+  4 1 — - — J  [  cos  2P  -  exp(-  p  tan  P)  cos(p  - 

2P)1 

(A-12) 

and  Qyjjgjjg  is  the  absorption  efficiency  from  van  de  Hulst^, 


exp(-  2p  tan  p)  exp(-  2p  tan  P)  -  1 
Q  =  1  +  - r -  + - : - 

vHabs  p  tan  P  ^ 


(A-13) 


2(p  tan  p) 


where 

and 

p  =  2X(n  -  l) 


(A-14) 


(A-15) 


A- 3 


In  Eqs .  A-10  and  A-11,  the  correction  term,  D,  is  given  by 


D  = 


[f(P)  +  1]  +  ,Zro-r  ^or  p^5(n-l) 


1.  632n 


|^[f(P)  +  l]P  for  5(n-l)^p^ 


(n-l)f(P) 


4.  08 


4.  08 


4.  08 


(A-16) 


(n  -  l)[f(P)  +  l]  .  _ - <  p  ^ - - — 

2n(l  +  3tanP)  1  +  3  tan  (3  1  +  tanP 


2.  04(n  -  l)[f(P)  +  1]  ^  _  4.08 

- - - -  for  P  >  - T 

nf(p)p  1  +  tanp 


where 


f(P)  =1+4  tan  p  +  3  tan^  p 


(A-17) 


At-4 
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ations  for  the  Ackerman  and  Stephens  Approximation 


The  Ackerman  and  Stephens  approximation  for  is  computed 


**'""*^  cos  psin(p  -  p)  +  — ^cos(p-2p) 


where 

and 

p  =  2X(n  -  l) 


(A-19) 


(A-20) 


(A-21) 


A- 5 


Equations  for  the  Mussenzveig  and  Wiscoiabe  Approximation^ 

The  Nussenzveig  and  Wiscombe  approximation  for 
computed  as 


2  -  1 


Q 


ext 


=  2  +  1.  9923861  X  %  8  ImjlCm^  +  l)  (m^  -  l)  X 

-  m'(m  +  1)' ‘  U' -  l)' ‘  [l  ^  ‘ 

-  1 

X  exp[2i(m  -  l)x]  -  7(m  -  l)  [^  “  (  2  )]  (  m  +  1 ) 


2j 


X  exp^2i  (m  —  1  +  2  jm)xj  |  0.  715353 7 X 


iff 


-  0.3320643  Im  Le  ^  l)  (m^+  l)(2m^  -  6m^  +3). 

+  o(x  +  ripple 

» 

Qabs  given  by 

where 


(A-22) 


(A-23) 


2  JL 


(Qabs)  =  I /Q^<t>O'lx)sin0  cosGde 

F  X  =  1 


-1  2  2 


(°.J_  =2  'x-3  1^4>(rj/)dx 


(A-24) 


(A-25) 


and 


i  2  2  X 


(0.J 

b,e,  ^ 


(A-26) 


A-6 


In  Eqs.  A-24  through  A-26, 


(A-27) 


where 


R 


3^ 


j,X,  =  l,2;  R^^=(-l)(2,-ue^)/( 


2  +  ue 


0  <A-28 


z  =  cos  6#  u  =  m  cos  0  ,  sin  0  =  m  sin  0 , 


(A-29) 


-2 

e^  =  l,  e^  =m  ,  2^  =  z, 


z  for  Eq.  A-24 

for  Eqs.  A-25  ‘.and  A-26 


(A-30) 


and  finally, 


b  =  4X  Im(m  cos  0  +  0’  sin  0) 


(A-31) 


Equations  for  the  Bohren  and  Nevitt  Approxiiation^^ 

The  Bohren  and  Nevitt  approximation  for  Q^j^g  is  computed  as 


4n' 


'  dbs 


(n  +  1)  -  (n  -  1)  exp(-t) 


J_  2_ 
2 

n  X 


^-xVn^-  1  ^ 

expV - n - J 


^  ^  1  -  exp(—  x)  (l  + 

n  / 


x) 


1 


(A-32) 


where 


X=  4kX 


(A-33) 
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Equations  for  the  Levine  Approximation^^ 

The  Levine  approximation  for  is  computed  as 


^abs  ~  -  A)lH(4ky)+  h(b)  -  H(4kX+  B)J 


{A-34) 


where  the  function  H(u)  is  given  by 


H(u)  =T->-  V--*- 


(A-35) 


and 


2n 


B  = 


1  + 


(n^-.kO 


1  -  A 


(A-36) 


The  value  for  A  in  Eqs.  A-34  and  A-36  is  determined  using  an 
iterative  scheme  given  by 


A  =  R  +  (a  -  l)  exp 


-  2nf 1 + 


n"  +  k 


1  -  A 


(A-37) 


where  the  iteration  is  begun  with  =  R  and 


R  = 


(n-  l)'  +  k' 

/  I*'  2 

( n  +  l)  +  k 


(A-38) 
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